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A combination of experimental techniques and molecular dynamics (MD) computer simulation 
is used to investigate the diffusion dynamics in AlgoNi2o melts. Experimentally, the self-diffusion 
coefficient of Ni is measured by the long-capillary (LC) method and by quasielastic neutron scat- 
tering. The LC method yields also the interdiffusion coefficient. Whereas the experiments were 
done in the normal liquid state, the simulations provided the determination of both self-diffusion 
and interdiffusion constants in the undercooled regime as well. The simulation results show good 
agreement with the experimental data. In the temperature range 3000 K> T > 715 K, the interdif- 
fusion coefficient is larger than the self-diffusion constants. Furthermore the simulation shows that 
this difference becomes larger in the undercooled regime. This result can be refered to a relatively 
strong temperature dependence of the thermodynamic factor <f>, which describes the thermodynamic 
driving force for interdiffusion. The simulations also indicate that the Darken equation is a good 
approximation, even in the undercooled regime. This implies that dynamic cross correlations play 
a minor role for the temperature range under consideration. 

PACS numbers: 64.70.Pf, 61.20.Ja, 66.30.Hs 
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I. INTRODUCTION 

Multicomponcnt liquids exhibit transport processes 
due to concentration fluctuations among the different 
components. In the hydrodynamic limit, these processes 
are described by interdiffusion coefficients [H,[2,[l]. In the 
simplest case of a binary AB mixture, there is one inter- 
diffusion coefficient -Dab- This quantity plays an impor- 
tant role in many phenomena seen in metallic mixtures, 
such as solidification processes [ij , the slowing down near 
the critical point of a liquid-liquid demixing transition Q 
or glassy dynamics [6|. 

Many attempts have been undertaken for different bi- 
nary systems to relate Dab to the self-diffusion con- 
stants -Da and Db via phenomenological formula (see 
e.g. 0, H3, HH, ) ■ An example is the Darken equa- 
tion 1 1 31 ] that is widely used to estimate the interdiffusion 
constant of simple binary fluid mixtures. This equation 
expresses -Dab as a simple linear combination of the self- 
diffusion coefficients, Dab = $(cbDa + caDb) (with c\, 
cb the mole fractions of A and B particles, respectively). 
Here, the so-called thermodynamic factor $ contains in- 
formation about static concentration fluctuations in the 
limit of long wavelength. 

The relationship between one-particle transport and 
collective transport properties is a fundamental question 
in undercooled liquids [6| . In the framework of the mode- 
coupling theory of the glass transition, Fuchs and Latz 
[l4j have studied a binary 50:50 mixture of soft-spheres 
with a size ratio of 1.2. Their numerical data indicate 
that the Darken equation is a good approximation for the 



latter system in the undercooled regime. However, from 
experiments or computer simulations, not much is known 
about the validity of the Darken equation for undercooled 
liquids. This is due to the lack of experimental data for 
interdiffusion coefficients in this case. Moreover, most of 
the computer simulation studies on the relation between 
self-diffusion and interdiffusion have been only devoted 
to the normal liquid state. In this case, the Darken equa- 
tion often seems to work quite well 0, E3, EH EE EH ■ 

In this work, a combination of experiment and molecu- 
lar dynamics (MD) simulation is used to study the diffu- 
sion dynamics in the metallic liquid AlsoN^o- In the MD 
simulation, the interactions between the atoms are mod- 
elled by an embedded atom potential proposed by Mishin 
et al. 11711 . The present work is a continuation of a recent 
study [18J, where a combination of quasielastic neutron 
scattering (QNS) and MD simulation was applied to in- 
vestigate chemical short-range order and self-diffusion in 
the system Al-Ni at different compositions. In the lat- 
ter study, we have shown that the MD simulation yields 
good agreement with the QNS data, both for structural 
quantities and the Ni self-diffusion constant, Dm- In the 
present work, an additional experimental method, the 
long-capillary (LC) technique, is used. This method al- 
lows to determine simultaneously the self-diffusion con- 
stant DNi and the interdiffusion coefficient Dab (see be- 
low). 

Above the liquidus temperature (i.e. in the normal liq- 
uid state) , thermodynamic properties as well as structure 
and dynamics of AlgoNi™ have been studied by different 
approaches (see, e.g., 0, El, HI HI, HI SI). The Al- 
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Ni system is an ordering system which is manifested in 
a negative enthalpy of mixing [24[. Thus, it does not 
exhibit a liquid-liquid miscibility gap where one would 
expect that the interdiffusion coefficient vanishes when 
approaching the critical point, whereas the self-diffusion 
constants are not affected by the critical slowing down 
(see (25| and references therein). Such a behavior is not 
expected for the system Al-Ni. 

In the computer simulation, the AlgoNi2o melt can be 
undercooled to an arbitrary extent avoiding the occur- 
rence of crystallization processes. Therefore, we were 
able to study a broad temperature range in our MD sim- 
ulations, ranging from the normal liquid state at high 
temperature to the undercooled liquid at low tempera- 
ture. In the experiments presented below crystallization 
occurs due to heterogeneous nucleation. Thus, the ex- 
periments were performed above the liquidus tempera- 
ture Tl ~ 1280 K. The combination of experiment and 
simulation presented in this work allows for a test of the 
validity of the Darken equation in AlgoN^o. We will see 
below that this equation is indeed a good approximation, 
even in the undercooled regime. 

In the next section, we summarize the basic theory on 
self-diffusion and interdiffusion. The details of the exper- 
iments and simulation are given in Sec. IHII and Sec. IIV1 
respectively. In Sec. [V] we present the results. Finally, 
we give a summary of the results in Sec. IVII 



Here, (t) is the position of particle j of species a at 
time t. Note that Eq. (|3|) is equivalent to the Green- 
Kubo formula ©. 

Interdiffusion is related to the collective transport of 
mass driven by concentration gradients. The transport 
coefficient that describes this process is the interdiffusion 
constant -Dab which can be also expressed by a Green- 
Kubo relation, i.e. by a time integral over an autocorre- 
lation function. The relevant variable in this case is the 
concentration or interdiffusion current [l[ given by 
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where ca = Na/N = 1 — cb is the total concentration 
(mole fraction) of A particles. As a matter of fact, the au- 
tocorrelation function of the variable Jab (t) depends on 
the reference frame and fluctuations of Jab (t) have to be 
adapted to the ensemble under consideration. Whereas 
experiments are usually done in the canonical ensemble, 
in a molecular dynamics simulation, the natural ensemble 
is the microcanonical ensemble with zero total momen- 
tum Hi. Thus, 
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II. SELF-DIFFUSION AND INTERDIFFUSION: 
BASIC THEORY 



follows, where toa and me denote the masses of A and 
B particles, respectively. Introducing the "centre of mass 
velocity of component a (a = A, B)" by 



Consider a three-dimensional, binary AB system of 
N = Na + Nb particles (with Na, Nb the number of 
A and B particles, respectively). The self-diffusion con- 
stant D SQ (a — A, B) is related to the random-walk mo- 
tion of a tagged particle of species a on hydrodynamic 
scales. It can be calculated from the velocity autocorre- 
lation function 
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via a Green-Kubo integral: 

/>oo 

D s . a = / C a (t)dt. (2) 
Jo 

In Eq. (fT]) , vj Q ' 1 (t) is the velocity of particle j of species 
a at time t. 

The self-diffusion constant can be also expressed 
by long-time limit of the mean-squared displacement 
(MSD): 
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we can use expression ([5]) to simplify the formula for the 
interdiffusion current, 

J A B(t)=iVCBC A fl + ^^ > ) V A (t). (7) 

V m B c B J 

Thus, we have to consider only the velocities of one 
species to compute Jab(0- 

Now, the autocorrelation function for the interdiffusion 
current is given by 



C AB (t) = (JabW-Jab(O)) 
= N 2 (c b ca) 2 (l + ' 
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The Green-Kubo formula for -Dab reads 
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where S cc (0) is the concentration-concentration struc- 
ture factor in the limit q — > 0. The function S cc (q) is the 
static correlation function associated with concentration 
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fluctuations. It can be expressed by a linear combination 
of partial static structure factors S a p(q) (a, (3 = A, B) as 
follows 0: 



S C c(q) 
with 
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Using elementary fluctuation theory [![, S cc (0) can be 
related to the second derivative of the molar Gibbs free 
energy g, 

c A c B d 2 g 
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S cc (q = 0) 

In Eq. [T2l &b is the Boltzmann constant and T the tem- 
perature. In the following, we will refer to $ as the ther- 
modynamic factor. 

We note that the total structure factor for the number 
density, S nn (q), and cross correlation between number 
density and concentration, S nc (q), can also be written as 
a linear combinations of partial structure factors. These 
functions are given by [l[ 

S nn (q) = S AA (q)+2S AB (q) + S BB (q) , (14) 
S„ c (q) = c B S AA (q) - c A S BB {q) + {c B - c A )S AB (q\lh) 

The typical behavior of these functions for a liquid mix- 
ture will be discussed in the result's section. The func- 
tions S nn (q), S nc (q) and S c Jq) are often called Bhatia- 
Thornton structure factors [271 ]. In principle, these func- 
tions can be determined in neutron scattering experi- 
ments, either by using isotopic enrichment techniques 
(see, e.g., Ref. |2l|) or by applying a combination of neu- 
tron scattering and X-ray diffraction [28| . 

With Eqs. (IS]) and ([13")) , the interdiffusion constant can 
be written as 
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Alternatively, -Dab can be also easily related to the self- 
diffusion constants to yield 
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where the functions A a p(t) denote distinct velocity cor- 
relation functions, 
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Note that the three functions A a p(t) can be expressed 
by the "centre-of-mass" correlation function Cab (t) and 
the velocity autocorrelation functions C a (t) (the latter, 
multiplied by l/c Q , has to be subtracted in the case of 
A AA (i) and Abb(*)) EH- Thus, the functions A a , 3 (t) 
do not contain any additional information compared to 
C AB (t) and C a (t) and so we do not consider them sepa- 
rately in the following. 

If one denotes the distinct part in (|T7|) by 
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one can rewrite Eq. (TT7|) . 



with 



D AB = <I>S(c a D b + c b D a ), 
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The quantity S measures the contribution of cross cor- 
relations to -Dab- If S = 1 holds, the interdiffusion con- 
stant is determined by a linear combination of the self- 
diffusion constants. In this case, Eq. ([2H)l leads to the 
Darken equation Note that, in the context of chem- 
ical diffusion in crystals, S is called Manning factor [29I ]. 

As in the case of self-diffusion, the interdiffusion con- 
stant can be also expressed via a mean-squared displace- 
ment which involves now the centre-of-mass coordinate 
of species A, 

, N A 

R A(*) = ^E r i A) (*)" ( 22 ) 
Then, the "Einstein relation" for Dab reads 

x [R A (i)-R A (0)r 

Dab = lim 1 + Nc A c B <f>- 

t^oo y m B c B J or 

(23) 

This formula can be used to determine Dab hi a com- 
puter simulation, where the system is located in a simu- 
lation box with periodic boundary conditions. However, 
in this case one has to be careful because the difference 
Ra(£) — Ra(0) has to be calculated in an origin inde- 
pendent representation [3(|. This can be achieved by 

computing this difference via the integral J Y A (t')dt' . 



III. EXPERIMENTAL METHODS 

A. Long— capillary technique 

The long-capillary technique (LC) has been used to 
measure interdiffusion and Ni self diffusion in liquid 
Al8oNi2o- The sample material production is similar to 
that of Alg7NiioCe3, which is described in Ref. [3l[. The 
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FIG. 1: Typical concentration profiles of a combined interdif- 
fusion and self-diffusion experiment. The squares denote the 
Al and Ni concentrations measured by energy-dispersive X- 
ray spectrometry (EDS) and the dots denote the 62 Ni concen- 
tration measured by inductively-coupled plasma mass spec- 
trometry (ICP-MS). The lines represent the best fit (least- 
square method) of the appropriate solution of Fick's diffusion 
equations to the measured concentrations. 



experimental apparatus, the measurement of the concen- 
tration profiles and the evaluation of the concentration 
profiles, including the determination of Fick's diffusion 
coefficients, are also described elsewhere [H, HI)]. Thus, 
here the experimental technique is reported only briefly. 
In more detail we describe an improved diffusion cou- 
ple setup, which has been used in this work. This setup, 
with a vertical diffusion capillary of 1.5 mm diameter, has 
an increased stabilization against natural convection and 
minimizes the systematic error of convective mass flow 
contributions to the total mass transport. 

The improvement of the diffusion couple setup implies 
the combination of intcrdiffusion and self-diffusion mea- 
surements in one experiment. An AI80M20 slice of 2 mm 
thickness, containing the enriched stable 62 Ni isotope, 
is placed between both rods of an interdiffusion couple. 
The interdiffusion couple consists of a 15 mm long rod of 
AI85N115, placed above the slice, and a 15 mm long rod of 
AI75M25, placed below the slice. This configuration al- 
lows the development of an error function shaped chem- 
ical intcrdiffusion profile simultaneously to the develop- 
ment of a Gauss function shaped self-diffusion profile. In 
a first approximation the diffusion of the enriched stable 
isotope takes place at the mean concentration AlgoN^o 
without influence of the changing chemical composition 
of the melt in the diffusion zone. The only necessary cor- 
rection results from the mass spectrometric measurement 
of the self-diffusion profile. Here the measured isotope 
incidences i( 62 Ni) of 62 Ni have to be corrected for the 
overlaying chemical concentration profile of natural Ni, 



cm, by using the following formula: 

c( 62 Ni) = c Ni (z( 62 Ni) - t( 62 m )) (24) 

with z( 62 Ni°) the natural incidence of 62 Ni and c( 62 Ni) 
the concentration of this Ni isotope with respect to all 
Ni isotopes. Typical concentration profiles of a diffusion 
experiment are given in Fig. [1] 

The diffusion couple configuration minimizes the risk 
of convection compared to conventional self-diffusion ex- 
periments in pure melts because of the solutal stabilized 
density profile of the melt column. This stabilizing ef- 
fect has been described in Refs. [34], [35[. In a standard 
self-diffusion experiment without chemical gradient the 
solutal stabilization effect is only due to the enrichment 
of a tracer. 

As a test for other mass transport processes wc mea- 
sured the mean-square penetration depth x 2 of interdif- 
fusion as a function of time t. We found a deviation from 
the linear behavior x 2 = 2£>ab£- This has been iden- 
tified as sedimentation of AI3M2 during solidification of 
the diffusion sample. This additional mass transport was 
simply corrected by subtracting this contribution as an 
off-set of the measured total mass transport. This pro- 
cedure adds a 5-10% error to the uncertainty of the dif- 
fusion coefficient. The total error in the long-capillary 
measurements of the self- and interdiffusion coefficients 
is about 30-40%. 



B. Neutron scattering experiments 

The second experimental technique used in this work is 
quasielastic neutron scattering. In this case, the AlgoNi2o 
alloy was prepared by arc melting of pure elements un- 
der a purified Argon atmosphere. The measurements 
were done at the time-of-flight spectrometer IN6 of the 
Institut Laue-Langevin. The standard Nb resistor high 
temperature vacuum furnace of the ILL exhibits a tem- 
perature gradient over the entire sample at 1800 K that 
was less than five degrees and a temperature stability 
within one degree. For the scattering experiment wc 
used a thin-walled AI2O3 container that provides a hol- 
low cylindrical sample geometry of 22 mm in diameter 
and a sample wall thickness of 1.2 mm. 

An incident neutron wavelength of A = 5.1 A yielded 
an energy resolution of 5E ~ 92 ficV (FWHM) and an 
accessible wave number range at zero energy transfer 
of g = 0.4 — 2.0A . Measurements were performed 
at 1350 K, 1525 K, 1670 K and 1795 K in 2 hour runs 
each. A run at room temperature provided the instru- 
mental energy resolution function. The scattering law 
S(q, 10) was obtained by normalization to a vanadium 
standard, accompanied by a correction for self absorp- 
tion and container scattering, and interpolation to con- 
stant wave numbers q. Further, S(q, u>) was symmetrized 
with respect to the energy transfer fkv by means of the 
detailed balance factor. 
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FIG. 2: Normalized scattering law of liquid Ai8oNi2o- The 
data at 300 K represent the instrumental energy resolution 
function. The line is a fit with a Lorentzian function that is 
convoluted with the instrumental energy resolution function. 
Diffusive motion of the atoms leads to a broadening of the 
quasielastic signal from which the Ni self-diffusivity can be 
obtained on an absolute scale. 



Figure [2] displays S(q,u>) at q = 1.0 A of liquid 
AlsoNi2o at 1525 K and the crystalline alloy at 300 K at 
q = 0.8 A . Diffusive motion in the liquid leads to a 
broadening of the quasielastic signal. The data were fit- 
ted with an Lorentzian function that is convoluted with 
the instrumental energy resolution function. From the 
full width at half maximum of the quasielastic line r a q- 
dependent diffusion coefficient D(q) can be computed via 
D(q) = T/(2hq 2 ). Towards small q incoherent scattering 
on the Ni atoms dominates the signal and the diffusion 
coefficient D(q) becomes constant yielding an estimate 
of -D s .Ni- Thus, the self-diffusion constant -D s ,Ni can be 
determined on an absolute scale [HI, H3] ■ 



IV. DETAILS OF THE SIMULATION 

For the computer simulations of the binary system 
AlgoNi2o, we used a potential of the embedded atom type 
that was recently derived by Mishin et al. [13] ■ In a re- 
cent publication [l8j ]. we have shown that this potential 
reproduces very well structural properties and the self- 
diffusion constant of Al-Ni melts at various compositions. 
The present simulations are performed in a similar way 
as the ones in the latter work: Systems of N = 1500 
particles (Nm = 300, A^ai = 1200) are put in a cubic 
simulation box with periodic boundary conditions. First, 
standard Monte-Carlo (MC) simulations in the NpT en- 
semble [HI were used to fully equilibrate the systems at 
zero pressure and to generate independent configurations 



for MD simulations in the microcanonical ensemble. In 
the latter case, Newton's equations of motion were in- 
tegrated with the velocity Verlet algorithm using a time 
step of 1.0 fs at temperatures T > 1500 K and 2.5 fs at 
lower temperatures. 

The masses were set to 26.981539 amu and 58.69 amu 
for aluminum and nickel, respectively. At each tempera- 
ture investigated, we made sure that the duration of the 
equilibration runs exceeded the typical relaxation times 
of the system. The temperatures considered were 4490 K, 
2994K, 2260K, 1996K, 1750K, 1496K, 1300K, HOOK, 
998 K, 940 K, 893 K, 847 K, 810 K, 777 K, 754 K, 735 K, 
715 K, 700 K, 680 K, and 665 K. In order to improve the 
statistics of the results we averaged at each temperature 
over eight independent runs. At the lowest tempera- 
ture, the duration of the microcanonical production runs 
were 40 million time steps, thus yielding a total simu- 
lation time of about 120 ns. The latter production runs 
were used to study the tagged particle dynamics. For 
the calculation of the intcrdiffusion constant Dab addi- 
tional production runs were performed in the tempera- 
ture range 4490 K > T > 715 K that extended the pro- 
duction runs for the tagged particle dynamics by about a 
factor of ten. This was necessary in order to obtain a rea- 
sonable statistics for Dab- Note that Dab is a collective 
quantity that does not exhibit the self-averaging prop- 
erties of the self-diffusion constant and thus it is quite 
demanding to determine transport coefficients such as 
the interdiffusion constant or the shear viscosity from a 
MD simulation. 



V. RESULTS 

In Eq. (j2"0|) , the interdiffusion constant Dab is ex- 
pressed as a linear combination of the self-diffusion con- 
stants. The prefactor in this formula is a product of 
the thermodynamic factor <E> and the Manning factor S. 
Whereas $ can be computed from structural input, the 
Manning factor contains the collective dynamic correla- 
tions in the expression for Dab (see Sec. [II]). In the 
following, we compare the simulated diffusion constants 
for AlgoNi2o to those from experiments. Moreover, the 
simulations are used to disentangle differences between 
self-diffusion constants and the interdiffusion constants 
with respect to the thermodynamic quantity $ and the 
dynamic quantity S. 

First, we discuss static structure factors at different 
temperatures, as obtained from the MD simulation. Fig- 
ure [5] displays the different partial structure factors at 
the temperatures T = 2000 K and T = 750 K. At both 
temperatures, a broad prepeak around the wavenumbcr 
q = 1.8 A -1 emerges in the NiNi correlations, which 
indicates the presence of chemical short-ranged order 
(CSRO). This feature is absent in the A1A1 correlations. 
In a recent work [l8T ] , we have found that the prepeak in 
<SNiNi(<?) is present in a broad variety of Al-Ni compo- 
sitions, ranging from xjsh = 0.1 to x-^i = 0.9. However, 
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FIG. 3: Partial structure factors, as obtained from the MD 
simulation, for a) T = 2000 K and b) T = 750 K. The multipli- 
cation by l/(c a C/3) 1//2 is introduced to increase the amplitude 
of SNiNi(?) relative to that of SaiaiC?). Note that the factor 
l/(c a c ( 3) 1//2 leads also to the asymptotic value S aa {q) = 1 for 
q — > co. 



the width of the prepeak decreases significantly with in- 
creasing Ni concentration and, in melts with a high Ni 
concentration, it appears also in Saiai(<?)- The prepeak 
in S a p(q) describes repeating structural units involving 
next-nearest a/3 neighbors which are built in inhomogc- 
neously into the structure. Of course, for the Al rich sys- 
tem AI80N120 considered in this work, only next-nearest 
Ni-Ni units exhibit the CSRO that is reflected in the 
prepeak. 

From the partial static structure factors, the Bhatia- 
Thornton structure factors can be determined according 
to Eqs. (p~0|) . (JT3J) and (Tl"5|) . These quantities are shown in 
Fig. gl again at T = 2000 K and at T = 750 K. Although 
these structure factors look very different for q > 2 A -1 , 
they are essentially identical in the limit q —> 0. As 
we have indicated before, the static susceptibility, asso- 
ciated with concentration fluctuations, can be extracted 
from the structure factor S cc (q) in the limit q — ► 0. As 
we can infer from Fig. 0] at the temperature T = 750 K 
the value of this susceptibility is very small. The small 
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FIG. 4: Bhatia-Thornton structure factors, as obtained from 
the MD simulation, at a) T = 2000 K and b) T = 750 K. 



value of S cc (q = 0) reveals that concentration fluctua- 
tions on large length scales are strongly suppressed. This 
is the typical behavior of a dense fluid that exhibits a 
strong ordering tendency. In contrast, at a critical point 
of a demixing transition a divergence of S cc (q = 0) is 
expected. 

As we have seen in Sec. [Til the ratio Dab/& can be 
expressed as a linear combination of the self-diffusion 
constants, provided 5 = 1 holds. In order to quantify 
the temperature dependence of S, we first define the fol- 
lowing mean-squared displacements: 



(r 2 (t)) int 



1 . m A c A \ at 

1 H Nc A c B x 

m B c B / 

([R A (t)-R A (0)] 2 ) 



(25) 



(r 2 (t)) scU = c A 



^-E(H B) W-rf(0)] 2 ) + 



r5 A) (0) ) (26) 



Whereas the intcrdiffusion constant can be calculated 
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FIG. 5: Simulation results of mean squared displacements 
(MSD) for self-diffusion (dashed lines) and interdiffusion 
(solid lines) for the temperatures T = 3000 K, 2000 K, 1500 K, 
1000 K, 850 K, 750 K, 715 K, and 665 K (corresponding to the 
curves from left to right. Note that for T — 665 K only 
(r 2 (i)) sc if was calculated. For the definitions of the MSD's 
see Eqs. (J25J and pBl. 



via D AB = lim t _ +00 $(r 2 (t)) int /(6t), the equation D AB = 
lim t _, 00 <i>(r 2 (i)) sc if/(6i) is only correct for 5 = 1. Figure 
[5]shows the quantities (r 2 (t))i nt and (r 2 (t)) se if for the dif- 
ferent temperatures. Both MSD's show a very similar be- 
havior. At high temperature, a crossover from a ballistic 
regime (oc t 2 ) at short times to a diffusive regime (oc t) at 
long times can be seen. At low temperature, a plateau- 
like region develops at intermediate times, i.e. between 
the ballistic and the diffusive regime. With decreasing 
temperature, the plateau becomes more pronounced. In 
(r 2 (i)) se if, the plateau indicates the so-called cage ef- 
fect p. The tagged particle is trapped by its neighbors 
on a time scale that increases with decreasing tempera- 
ture. Although the MSD for the interdiffusion, (r 2 (t))i nt , 
describes also collective particle transport, the plateau 
in this quantity has the same origin: The particles are 
"arrested" on intermediate time scales. Moreover, the 
differences between (r 2 (t)) sc \f and (r 2 (t))i nt are anyway 
very small in the whole time and temperature range un- 
der consideration. This means that the cross correlations 
do not give a large contribution to (r 2 (i))i n t. 

From the MSD's in Fig. [51 the Manning factor S can 
be extracted using Eq. (|2Tj) . In Fig. [6] we see that the 
Manning factor varies only slightly over the whole tem- 
perature range, located around values between 0.8 and 
1.0. Also shown in Fig. [fj] is the thermodynamic factor i> 
and the product $5. We have extracted $ from the ex- 
trapolation of the structure factors S cc (q) toward q — > 
[see Eq. (|13p]. In contrast to the Manning factor 5, the 
thermodynamic factor $ increases significantly with de- 



creasing temperature and thus, also the change in the 
product $5 is dominated by the change in <f>. There- 
fore, differences in the qualitative behavior between the 
self-diffusion constants and the interdiffusion constant 
are dominated by the thermodynamic factor. 

An Arrhenius plot of the diffusion constants as ob- 
tained from simulation and experiment is shown in Fig. [7] 
The self-diffusion constants Z?Ni and -Dai from the simu- 
lation are very similar over the whole temperature range 
4490 K> T > 665 K. In a recent publication [3, we have 
found that, in the framework of our simulation model, 
this similarity of the self-diffusion constants occurs in Al 
rich compositions of the system Al-Ni, say for ca\ > 0.7. 
Whether this is also true in real systems is an open ques- 
tion. However, the neutron scattering results for D^i as 
well as the single point obtained from the LC measure- 
ment is in very good agreement with the simulation data. 



Asta et al. l_5j have computed the concentration de- 
pendence of the self-diffusion constants at T = 1900 K 
using two different embedded atom potentials, namely 
the one proposed by Voter and Chen [3^] and the one 
proposed by Foiles and Daw [4(j. For both potentials, 
they find very similar values for L>Ni and Dai in AlgoNi2o, 
in agreement with our results. However, their results for 
the Ni diffusion constant are significantly higher than 
the ones found in our quasielastic neutron scattering ex- 
periment and our simulation. They report the values 
L> Ni w 1.5 ■ 10" 8 m 2 /s and D Ni w 1.9 • 10~ 8 m 2 /s for the 
Voter-Chen potential and the Foiles-Daw potential, re- 
spectively, whereas we obtain Dni ~ 10~ 8 m 2 /s from sim- 
ulation and experiment. Thus, the potential proposed by 
Mishin et al. [l7j|, which is used in this work, leads to 
a better agreement with the experiment, as far as self- 
diffusion in Al8oNi2o is concerned. 

We emphasize that the statistical error in both the 
neutron scattering data and the simulation data for the 
self-diffusion constants is relatively small. In both cases, 
the error bars for the corresponding data points in Fig.0 
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FIG. 7: Arrhenius plot of interdiffusion and self-diffusion con- 
stants, as obtained from experiment and simulation, as indi- 
cated. The experimental results are measured by quasielastic 
neutron scattering (QNS) and by the LC technique. The lines 
through the data points are guides to the eye. The vertical 
dotted line in a) marks the location of the experimental liq- 
uidus temperature, TL ~ 1280 K. The vertical dashed line is 
at the location of the critical temperature of mode coupling 
theory, T c « 700 K, as estimated by the MD simulation [4lj . 
Panel b) is an enlargement of the data of panel a) in a tem- 
perature range above Tl. The error bars of simulation and 
QNS data are of the order of the size of the symbols. 



are smaller than the size of the symbols. 

Due to the lack of self-averaging, it is much more dif- 
ficult to yield accurate results for Dab from the simu- 
lation. Therefore, in this case we considered a smaller 
temperature range than for the self-diffusion constants 
to yield results with reasonable accuracy. As we can in- 
fer from Fig. [71 the interdiffusion constant is larger than 
the self-diffusion constants over the whole temperature 
range. The difference becomes more pronounced with 
decreasing temperature. At T = 715 K, the diffusion co- 
efficient Dab is about a factor of 3 larger than D^i and 



Dai- This behavior is of course due to the increase of 
the thermodynamic factor $ at low temperature. Also 
included in Fig. [JJ are the results of the LC measurements 
of .Dab and DNi- These results are much less accurate 
than those of the quasielastic neutron scattering experi- 
ments for the determination of Z?Ni (see the error bars for 
the LC data in Fig.UJ)). Nevertheless, the LC data show 
that -Dab > ^Ni holds, in agreement with the simulation 
results. 



VI. CONCLUSION 

A combination of experiment and molecular dynamics 
(MD) simulation has been used to investigate the diffu- 
sion dynamics in liquid AlsoNi2o- We find good agree- 
ment between simulation and experiment. Both in ex- 
periment and in simulation, the interdiffusion constant 
is higher than the self-diffusion constants. This is valid 
in the whole temperature range considered in this work, 
i.e. in the normal liquid state as well as in the undercooled 
regime. In the latter regime (which is only accessible by 
the simulation), the difference between the interdiffusion 
constant and the self-diffusion constants increases with 
decreasing temperature. 

All these observations can be clarified by the detailed 
information provided by the MD simulation. Both the 
thermodynamic factor $ and the Manning factor S have 
been estimated directly and accurately over a wide tem- 
perature range, as well as self-diffusion and interdiffusion 
coefficients. The central result of this work is shown in 
Fig. [51 Whereas the thermodynamic factor $ increases 
significantly by lowering the temperature, the Manning 
factor S shows only a weak temperature dependence. 
Moreover, the value of S is close to one which means that 
dynamic cross correlations are almost negligible and thus, 
even in the undercooled regime, the Darken equation is 
a good approximation. The temperature dependence of 
$ is plausible for a dense binary mixture with a strong 
ordering tendency. The situation is similar to the case of 
the isothermal compressibility which normally decreases 
with temperature in a densely packed liquid leading to 
very low values in the undercooled regime. In the same 
sense, the response to a macroscopic concentration fluc- 
tuation described by S cc (q = 0) tends to become smaller 
and smaller towards the undercooled regime which corre- 
sponds to an increase of <E> with decreasing temperature 
(since $ oc 1/S cc (q = 0)). 

We note that the data shown for Dab are all above the 
critical temperature T c of mode coupling theory which 
is around 700 K for our simulation model (see Fig. [7]) 
[4lj . Since it is expected that the transport mechanism 
changes below T c [a], it would be interesting to see how 
such a change in the transport mechanism is reflected in 
the interdiffusion constants. This issue is the subject of 
forthcoming studies. 
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